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Development of a Computer Algorithm for the Analysis of Variable-Frequency 

A.C. Drives - Case Studies Included 


M. David Kankam 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

and 


Owen Benjamin* 
Cornell University 
Ithaca, New York 14850 


ABSTRACT 

This report documents the development of a computer software for performance prediction 
and analysis of voltage-fed, variable-frequency AC drives for space power applications. The AC 
drives discussed include the pulse width modulated inverter (PWMI), the six-step inverter and the 
pulse density modulated inverter (PDMI), each individually connected to a wound-rotor induction 
motor. Various d-q transformation models of the induction motor are incorporated for user- 
selection of the most applicable model for the intended purpose. 

Simulation results of selected AC drives correlate satisfactorily with published results. 
Future additions to the algorithm are indicated in the Summary. These improvements should 
enhance the applicability of the computer program to the design and analysis of space power 
systems. 

L INTRODUCTION 

Variable-frequency drives have been extensively documented in the literature. Most 
publications deal with a particular problem area and application. In the realm of space mission 
applications, there is a need to incorporate flexibility for options and assessment on various 
potential voltage-fed inverters as possible drive sources for selected A.C. machinery. 

This report describes the development of a computer algorithm for performance prediction 
and the analysis of variable-frequency A.C. drives. These drives have a potential application in 
space missions. Among the AC drives discussed are the pulsewidth modulated inverter(PWMI), 
the six-step inverter and the pulse density modulated inverter(PDMI). A general synthesis of the 
PWMI output is described. The outputs of the inverters are d-q transformed for subsequent 
excitation of a similarly transformed, symmetrical induction machine model. The algorithm 
accommodates transformation in the arbitrary, fixed, synchronously rotating and multiple reference 
frames of operation. A complete flow-chart of the algorithm is provided in the Appendix. 

Simulation results obtained by the algorithm are compared with pertinent publication for 
verification. Salient characteristics of the results are discussed. Included in the summary are 
planned additions for subsequent incorporation in the algorithm. These extensions should further 
improve the usefulness of the algorithm for studies relating to space power applications. 


*1991 Summer OA1ZNASA Collaborative Aerospace Intern at Lewis Research Center. 



LL MODELING OF DRIVE SOURCES 


This section discusses various AC drives individually applied to a three-phase symmetrical 
induction machine. These drives permit a comparison of the effects of their output harmonic 
content on the performance characteristics of the induction machine. Also, the use of these drives 
demonstrates the flexibility of the developed algorithm in its choice of a particular drive for the 
intended application. 

ULL Pulse Width Modulated Inverter IPWMIl 

A Pulse Width Modulated Inverter has been simulated from development of an analytical 
procedure for deriving input harmonic voltages applicable to a variable frequency AC drive. The 
synthesis for obtaining an analytical expression for the line-to-line synchronized voltage 
waveform, v L . L (t), from a line-to-line voltage envelope, v e (t), and a continuous pulse 

operator, Vp(t), is shown schematically in Fig. II- 1. A summary of the method follows. The 
complete modeling details of the PWMI is contained in a NASA-Lewis 1990 Summer Report.fl] 



Fig. II- 1 - Synthesis of Line-to-Line Synchronized Voltage Waveform 


The voltage envelope controls the output frequency of the synchronized waveform, while 
the voltage pulse operator controls the number of pulses within a half cycle and the duration of time 
the switches stay closed. 

Reference [2] shows that the general synchronized voltage waveform, in Fig. II- 1, has the 
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following constraints: 


V > 6 

* (I) 

0<t/ t < 1 

P 

«/ dt = o 

where, 

E - amplitude (d.c. volts) of the line-line voltage envelope v e (t); 
u - 2 tt/ t * fundamental frequency of the output voltage (rad/s) 
of the inverter; 

w p “ 2 tt/T p » multi -pulse frequency (rad/s) 

* time interval between leading edges of two consecutive pulses 
of the voltage operator v p (t). 
r « pulse width (secs) of v p (t). 

T m - mode period (secs) - T/ 6 . 


The synthesized line-to-line output waveform of the PWMI is obtained from the relation: 
V L-LW = v e(*) * v p(0 (2) 


The voltage operator is defined in eqn. (3), 


r 


V -\ 


Ep ( t) 


0 < t < 7T /, 


3d) 

w/ 3«. < * « 2 ’' 3u 


gp(t-n/ 3|ij ) , 


(3) 


’P 3u>' ’ ‘ "'3w a) 

The Fourier expression for the unity amplitude multipulse operator is stated in eqn. (4). 
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( 4 ) 


Equation (5) shows the general Fourier expression of the line-to-line voltage 
V L-L (t) = l a k Coskut + b k Sinkwt, k=l,5,7,.. 


where the coefficients are 
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* | .Cosn+ZkSlng.CGjjj 
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Other variables in eqn. (6) are given in eqn. (7). 
F n = H* •'i-t-'osnu.pT 


H n =(^E) - k 2 

' CO ' 

G n = -2Sin/^ + + Sinn 

' 3to / 


J = -2Cos 
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/mrw , \ . „ 

( P + n) + Cosn 
'3 Li / 


( 5 ) 
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( 7 ) 


k = 1,5,7,11,13 
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Typical values used in subsequent simulation are summarized below. 
E = 115 v D.C 

a) = operating frequency 

= 12 Hz, 30 Hz, 54 Hz 

w eb = k ase frequency (rated) 

= 60 Hz 

io = 36 0 Hz 

P 



V, = line-to-neutral r.m.s voltage at u>=u> 
base ' 

= 51.77v 


11.2. Six-Step Inverter 


Figure II-2a shows a schematic for the basic circuit of a six-step inverter induction motor 
drive. Typically, the control of the magnitude of the dc link voltage, Vj, is achieved by 

maintaining a constant v/f ratio. 
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Fig. II-2a - Inverter-Induction Motor Drive 


The six- step inverter is designed so that the output separates in six distinct 60 degree periods in one 
cycle. This form of output is sometimes referred to as the “staircase” waveform, depicted in 
Fig. II-2b for one phase voltage of the inverter.[3] 
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Fig. II-2b - Output Voltage Waveform of a Six-Step Inverter 
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Fig. II-2b - Output Voltage Waveform of a Six-Step Inverter 


IL1 Pulse Density Modul ated Inverter (PDMU 

Analysis of the output of the Pulse Density Modulated Inverter is based on a synthesized 
waveform in reference^] and shown in Fig. II-3. Using a constraint of zero voltage switching, a 
low frequency voltage component, vpf, is derived from complete half cycles of the high - 
frequency link voltage, vp^. 
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Fig. II-3 - PDM Synthesis of low-frequency Voltage 

The third waveform, vpp, displayed in Fig. II-3 is the output synthesized from the first 

two input wave signals. This waveform may be Fourier analyzed and applied to a connected AC 
machine. This is achieved by either a programmed logic or some form of a modulation scheme. 
Founer analysis of the this particular frequency output is shown below. 

v LfO) = ^f a k * si n(k*ax) + b^ * cos(£*etx)}, 

where 
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a k = ( 2 * E m* 0) Lf) /71 * 
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and 

f c (t) = sin(tOHf*t)*cos(k*ci)Lf*t), 
d = 7i/(19*to L f). 

Expansion and simplification of eqn. (8) results in the expression of eqn. (9). 

a k = * { [ (-cos(d*a) - cos(3*d*a) + cos(4*d*a) +cos(6*d*oc) - cos(7*d*a) - 

cos(l l*d*a) + cos(12*d*a) + cos(14*d*a) - cos(15*d*a) - 
cos(17*d*a) + cos(18*d*a) - cos(19*d*a) + 2*cos(2*d*a) - 
2*cos(5*d*a) + 2*cos(8*d*a) - 2*cos(9*d*a) + 

2*cos(10*d*a) - 2*cos(13*d*a) + 2*cos(16*d*a))/a ] + 

[(-cos(d*P) - cos(3*d*p) + cos(4*d*P) + cos(6*d*p) - cos(7*d*P) - 
cos(l l*d*(3) + cos(12*d*P) + cos(14*d*P) - cos(15*d*P) - 
cos(17*d*p) + cos(18*d*p) - cos(19*d*P) + 2*cos(2*d*p) - 
2*cos(5*d*P) + 2*cos(8*d*P) - 2*cos(9*d*p) + 

2*cos(10*d*P) - 2*cos(13*d*P) + 2*cos(16*d*P))/ p ] }, (9) 

where 

a = toptf + (k * a>L f ) 

P = topjf - (k * a>Lf) 

A similar derivation holds for if 
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f s (t) = sin(t0Hf*t)*sin(k*C0Lf*t) 

is substituted for f c (t) in eqn. (8). This yields the expression of eqn. (11). 


( 10 ) 


bfc = [E m *G>Lf/rc] * {[( _ sin(d*a) - sin(3*d*a) + sin(4*d*a) + sin(6*d*a) - sin(7*d*a) - 
sin(l l*d*a) + sin(12*d*a) + sin(14*d*a) - sin(15*d*a) - 
sin(17*d*a) + sin(18*d*a) - sin(19*d*a) + 2*sin(2*d*a) - 
2*sin(5*d*a) + 2*sin(8*d*a) - 2*sin(9*d*a) + 

2*sin(10*d*a) - 2*sin(13*d*a) + 2*sin(16*d*a))/a ] + 

[(sin(d*P) + sin(3*d*P) - sin(4*d*P) - sin(6*d*P) + sin(7*d*P) + 
sin(ll*d*P) - sin(12*d*P) - sin(14*d*P) + sin(15*d*p) + 
sin(17*d*P) - sin(18*d*P) + sin(19*d*p) - 2*sin(2*d*p) + 

2*sin(5*d*p) - 2*sin(8*d*p) + 2*sin(9*d*p) - 

2*sin(10*d*p) + 2*sin(13*d*P) - 2*sin(16*d*p)) / p ] ) (11) 

Each of the above voltage sources is subsequently used as a driving source to simulate the 
performance of a wound-rotor induction motor. 

UL INDUCTION MACHINE MODELING 

This section describes the derivation of the d-q transformation model of an induction 
machine. Various forms of the model are discussed for the arbitrary, synchronous rotating and the 
multiple reference frames. 

IIL1 Basic Machine Equations 

A schematic diagram of the connections and conventions of an ideal, two-pole, three- 
phase, symmetrical induction machine and its magnetic axes representations are shown in Figs. III- 
1 and III-2, respectively. 



Fig. III-l - Connections and Current Conventions 

for Three-Phase Induction Motor 
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Fig. III-2 - Relative Positions of the Magnetic Axes of 

Stator and Rotor Windings. 


With the appropriate subscript for the stator (as, bs, cs) and rotor (ar, br, cr) phases, the following 
voltage eqn. (12) is applicable to each of the six windings shown in Fig.III-2. 

v = p* X + r* i (12) 

where X, r and p are the total flux linkage, the winding resistance and the differential operator dldt, 
respectively. 

The flux-linkage equations[4] for a three-wire machine are summarized in eqn. (13) 


L v CO.^r 


" co { e ' + H) 

L " co s (*'- 2 f) 


L„ cos 


V L* u 

L„ cos^e, - ^ L„ cos^fi, + 

L„ cosJ $ 7 Ljr cos M) 

L„ cos^e, + cose, 

L„ cose, L„ cos^e, + L„ cos^fl, 

Isfr COS^e, Lit COSe, Lgr COS^e, -f* — ^ 

L„ cos^e, + y'J cos^e, - j') L„ cose, 


L ss ~^s~ L sm 


L/r ~L r - L rm . 
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The parameters L sm and L rm are the mutual inductances between the individual stator and 

individual rotor phases, respectively. The term 6 r is the angular displacement between 

corresponding phases of the stator and rotor axes, as shown in Fig. ID-2. The time varying 
coefficients in eqn. (13) may be eliminated by means of transformation equations. 

11L2. Transform ation to d-q Variables 

Fig. IH-3 illustrates the angular relation of the stator and rotor axes relative to an orthogonal 
set of d-q axes. The transformation axes are rotating at an arbitrary electrical angular velocity, to. 

It is noteworthy that the stator as-bs-cs set is fixed in the stator. Similarly, the ar-br-cr set is fixed 
in the rotor and, thus, rotates at the rotor electrical angular velocity, o^.. 



Fig. III-3 - Axis of 2-pole, 3-phase symmetrical machine. 


In conjunction with Fig. III-3, and by means of geometric projections, it can be shown that 
the transformation from the machine a, b, c co-ordinates to the d-q reference frame is given by 
eqn. (14), 


fqs 
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(14) 
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where P = 0 - 0 r ( 15 ) 

and 0 is a continuous function. In these equations, the variable *f* represents either voltage, 
current, or flux-linkage. 


III.3 Reference Frames 

III. 3.1 Arbitrary Reference Frames 

The transformation equation (14) may be used to transform the stator and rotor voltages and 
currents of eqns (12) and (13) to the form shown in eqn. (16). 


v ds 


r s +L ss *p 

-o/ f L ss 

v qs 


of f L ss 

r s + ^ss*P 

( 

v dr 


M*p 

-(d)- (0 r )M 

1 


(co- to r )M 

M*p 


M*p 

- (t)*M 


l ds 


M*p 

* k 

i 

l qs 

/ / 

/ 


./ 

r r \L n *p 

-((O-wJL". 


l dr 

(o> co fjLfj. 

r r +L 'rr*P 


i' 

l qr 


(16) 


By setting the speed to, of the arbitrary reference frame, in eqn. (16), equal to the speed of 
the desired reference frame, the voltage equations may be expressed in any reference frame. The 
stationary or stator reference frame voltage equations may be obtained by setting co equal to zero in 
eqn. (16). Setting to equal to (Op the electrical angular velocity of the rotor, gives a reference frame 
fixed in the rotor. 

III. 3. 2 Synchronous Rotating Reference Frame 

The synchronously rotating reference frame voltage equations are obtained by setting (0 
equal to the electrical angular velocity, (0 e , of the fundamental frequency components of the 
applied stator voltages. This substitution yields equation (17). 
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v ds 


r S + ^SS*P 
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r s +L ss*P 
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All entries in the impedance matrix of either eqn. (17) or (16) are in actual units. For 
steady-state operation, the operator ‘p’ is replaced by the appropriate ‘ jco a ’ term where « a is the 
frequency of the voltages in the transformed frame. 


111.3.3 Multiple Reference Frame 

The voltage equations are linear if constant rotor speed is assumed. In that case, the 
principle of superposition can be used to obtain the machine performance due to the total harmonic 
content in the source voltages. For balanced sinusoidal applied voltages, the synchronously 
rotating d-q transformed voltages become constant voltages. This can be verified via application of 

eqn. (14) to such input source voltages, with the arbitrary speed o> set to the fundamental source 
frequency 0 ) e . The resulting constant voltages facilitate computation of constant- speed electrical 
transients and steady state performance by means of d.c. circuit theory. 

In the case of possible harmonic content in the source voltages, each of the harmonic 
components forms a balanced set of voltages in the synchronously rotating reference frame. E ac h 
of these balanced set of voltages becomes constant when transformed to a particular reference 
frame in which d.c. circuit theory can be used. If this is done, the resulting currents can then be 
back-transformed to the synchronously rotating reference frame. Subsequently, the final currents 
can be used to generate the constant-speed performance of the connected machine. Thus, the so- 
called multiple reference frame [4] application expedites the calculation of steady-state performance 
from any periodic stator voltage waveforms. Also, it eliminates the need to use phasors and 
complex impedances. 

The synchronously-rotated harmonic voltages, to be derived later, will show that each of the 
positively rotating voltages(that is fundamental, seventh harmonic, etc.) transforms to a voltage of 

frequency (k-l)co e , where k is the harmonic order. Thus, unity value of k denotes a d.c. or 
constant value for the transformed fundamental frequency component. Similarly, all negatively 
rotating voltages (k=5, 1 1, etc.) transform to voltages with frequency (k+l)o) e . Hence, the 
voltage equations for the symmetrical induction machine in any reference frame is given by eqn. 

(16) in which the arbitrary reference frame speed, 0), is set equal to the speed of the desired frame. 
When this is done, the resulting multiple reference frame equations are those summarized in eqn. 
(18). It is emphasized that in eqn. (18), the letter V takes on the value +k or -k, depending on the 
application of a positive or negative-sequence voltage, respectively. Since d.c. theory applies to all 
harmonic terms in the respective multiple reference frames, the time-derivative operator ‘p’ is set to 
zero in the steady-state performance calculation. 
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The impedance matrix, [Z], of eqn.( 18) has been used for computer simulation to verify 
some published results. For ease of comparison, the [Z] matrix has been per-unitized and 
expressed in eqn. (19). 




-fr*n*X ss 
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r s 

fr*n*X m 

/ 

0 


0 
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where the per-unit values in eqn. (19) are defined below, namely 
f r=0) e /( % 

fr* x ss = ( C[ V 0) b)* c,) b*^55 
f r * x m = 

fr*X / rr =(0) e /(%)*(%*L r r 
akf= n-((£> r / (0 e ) 

(n (O e - (O r )M = ( n - co r / 0) e )fr*X m = akf*fr *X m 
Mp =;'(A:-+l)*(ty e /a) b )*(O b *M = j(k-+\)fr*X m 


U1A Derivation of Multiple Reference Frame Voltages 

Consider a three-phase machine operated from a balanced source with a fundamental 
frequency, co e . Using eqn. (13), the d-q stator voltages fixed in the stator (0 = 0) are given by eqn 
( 20 ). 
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vS ds 
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i.e. 

,s . _ 


v ds = l/sqrt(3)(v bs - v cs ) 


v* = V 

y qs v as 


For a harmonic-rich source, the input voltages are periodic and sinusoidal. Hence, eqn. (21) gives 
the Fourier expansion of the transformed voltages. 

vS ds =|/ v kdc cos ^ ft) e r ) + v kds sin(/:cu e r)} 

v V = £/ v kqc cos (*' «V) + v kqs sin ^ w e r ) ) (21) 

When expressed in terms of the voltages v s ds and V s qs , the applied stator voltages in the 
synchronously rotating frame are shown in eqn. (22). 

v& ds = vS qs sin (°y) + v s d s cos(6) e t) = v +e kds + V e kds 

v V = v V C0s(6) ^) - v s ds sm(co e t) = v +e kqs + v' e kqs (22) 

Substitution of eqn. (21) into eqn. (22) will show that the voltages v +e kds and v +e kqs of the 
positive sequence applied voltages form a series of balanced sets. Fork>l, these sets rotate at an 
electrical angular velocity of (k-1 )co e in the counterclockwise direction with regard to the 

synchronously rotating frame. Similarly, v~ e kds and v e kqs of the negative sequence applied 

voltages rotate at the frequency -(k+1 )O e in the clockwise direction with respect to the 
synchronously rotating reference frame. 

In a manner consistent with the forms in eqn. (22), the “+-ke” variables in the multiple 
reference frame are obtained from the “+-e” of the synchronously rotating reference by use of 
eqns. (23) and (24). [4], 




COS {(A - \)(O e *t } 

sin{(A - 1 )co e *t }l 

' 



f +ke 
Jqs 

fdr +ke 


-sin{(jfc - l)(O e *t) 

cos{(* - 1 )(O e *t )l 
1 

[0] 

cos { (k - 1 )<o*t J sin{(* - 1 )(o e *t } 

* 

fkqs +e 

f kdr +C 

e +ke 

J qr J 



[0] 1 

-sin{(* - 1 )<o e *t) cos{(Jt - \)(o g *t } 


Jkqr * 


14 



7 1 

*• 

l 


COS {(it + 1 )(O e *t } 

-sin {(1 + l)<u *t }l 



~fkdi e ~ 

f -ke 
Jqs 


sin{(A + l)(O e *t) 

cos {(A + 1 )(O e *t }l 

[0] 

* 

fkqi 6 

4“ 



1 

cos { (1 + 1 )0) e *t } -sin{(A + 1 )(O e *t } 


f'kdr 6 

r / -ke 
jV 



[0] 1 

sin {(A + 1 )co*t) cos {(A + 1 )u>*t ) 

c j 

i 

ft ' e 

\Jkqr j 


( 24 ) 


+-e 


wher Qfd s + ' ke > fq S + ~ ke represent variables in the multiple reference frame, and/^ + ' e , /^ 5 " 

represent variables in the synchronous reference frame. It is noteworthy that, for stator only-fed 
induction machine, the rotor windings are shorted. Hence, for voltage representation, v ar, v br, 

and v cr are 0- Consequently, so are all the rotor variables subscripted “dr”, etc. in eqns. (23) and 
(24). 

Eqns. (25) represent a given set of 3 phase, balanced nonsinusoidal voltages. 


= £{ V £*cos (k* co e *t ) } , 
v bs = Z{V k *cos [k(<D*t - 2jc/3)] } , 

v cs = E{ V £*cos [k(co e *t + 27t/3)] } , 
where k = 1, 5, 7, 11,... 

The synchronously rotating reference frame voltages are as expressed in eqn. (26), 


(25) 


vC kds + ~ 


-+V/ c+ _ * sin((A -+ 1 )(O e t) 

v6 kqs + ~ 


V k+ _ * cos ((£ -+ l)(O e t ) 

* P -\ — 

v kdr 


0 

* P + - 

_ v kqr 


0 


(26) 


where the positive superscripts refer to positive-sequence applied voltages, k + = 1,7, 13... and 
the negative superscripts denote negative-sequence applied voltages, k. = 5, 11... The variable 
Vj c+ . represents the magnitude of the appropriate source harmonic voltage applied to the machine. 

Comparison of eqns. (23) and (24) and (26) shows that, for voltage considerations, eqn. (27) 
holds. 
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f +e kqr\ T =l' >e dsk+ v e qsk+ v ' 6 drk + v ' e qrk^ T 


^ +€ kds f +C kqs f+ £ kdr 

( 27 ) 

\f € kds f 6 kqs f € kdr ^kqr^ = ^dsk- v6 qsk- v ' € drk- v ' € qrk-^ 

Thus, substitution of eqn. (26) into eqns. (23) and (24) gives the multiple reference frame voltages 
in eqn. (28). b 


v ds + ' ke 


0 

V +-ke 

v qs 


+ 

1 

* +-kp 

v dr 


0 

v' +-ke 

_y 


0 


The impedance matrix of eqn. (17) and the voltage magnitudes of eqn. (26) may be used to 
simulate the machine performance in the synchronously rotating reference frame. 


For steady state operation, the operator ‘p’ in the impedance matrix is replaced by 

for stator applied positive sequence voltages, and “j(k_+l )(Qg” for negative sequence 

voltages. For simulation in the multiple reference frame, eqns. (19) and (28) are used. It may be 
recalled from section III.3.3 that all p-dependent terms are non-existent due to the resulting time- 
invariant voltages in the multiple reference frame. The following section discusses some 
simulation results. 


UL SIMULATION RF.STTI.TR 


Results of computer simulations of a symmetrical induction machine are described for the 
various drive sources discussed in section n. The developed computer algorithm for the 
simulation is flow-charted in Appendix Al. 

ULL Parameters of Induction Machine 


The following values in Table I are used for simulating the impedance matrix in the multiple 
reference frame. All values are in per-unit[4] unless otherwise stated. 
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TABLE I - Summary of Induction Machine Parameters 


Parameter 

Value* 

r s 

0.0453 

^ r 

0.0222 

x ls 

0.0775 

x 'lr 

0.0322 

X m 

2.0420 

x ss = X ls + x m 

2.1195 

X 'rr= X 'lr + X m 

2.0742 

(Ob 

2rc*60 rad/s 

a" 

II 

1.0 

s 

0.03 

(Of. = (0*{\ - s) 

(2tc* 60)(0.97) rad/s 


s - slip ; (Ofj - base electrical angular velocity 
i- per-unit values unless stated otherwise 


IV. 2. Six-Step Induction Motor Drive 

The simulation of a six-step induction motor drive seeks to validate the derived analytic 
relations and computer algorithm by comparing the simulated results with published data. The 
source input voltages are taken from reference[4]. The harmonic voltage magnitudes in the 
multiple reference frame are indicated in Table II. Note that the rotor voltages are zero due to 
shorting of the rotor windings. Also, consistent with eqn. (28), only the q-axis voltages exist in 
the multiple reference frame. 

TABLE II -Six-Step Voltage Magnitudes in Multiple Reference Frame 


Harmonic Order 


Value(pu) 


v 

V 

V 

V 

V 


+le 

qs 

-5e 

qs 

+7e 

qs 

-lie 

qs 

+13e 

qs 


4/tc 
4/5 7t 
-4/7 n 
- 4/ 11 n 
4/13 n 


The inverter constant input voltage V, set to 1 .0 pu for the simulation, is the amplitude of the 
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balanced squarewave voltages at the source. 

Table III shows a comparison between the published[4] and calculated steady-state currents 
in the multiple reference frame. The current values are corrected to the third decimal place. With 
the exception of the fundamental frequency components, the remaining calculated currents correlate 
reasonably with the published results. The discrepancy in the fundamental frequency currents is 
discussed later in conjunction with line current and electrical torque plots. 


TABLE III - Comparison of Calculated and Published Per-Unit d-q Currents 


n 

(+-k) 

,• ne 

Rgflil 

ds 

Calc. 

Rgfl4I 

: ne 

1 qs 

Calc 

«• ne dr 
RtfMl 

Calc. 

/ ne 

* qr 

RefI41 

Calc. 

+1 

0.600 

0.74120 

0.013 

1.50260 

0 

-0.1791 

0 

-1.5432 

-5 

-0.460 

-0.4602 

0.053 

0.053 

0.453 

0.453 

-0.053 

-0.0533 

+7 

-0.2360 

-0.2360 

-0.022 

- 0.022 

0.232 

0.232 

0.022 

0.022 

-11 

0.096 

0.096 

-0.005 

-0.005 

-0.095 

-0.095 

0.005 

0.005 

+13 

0.069 

0.069 

0.003 

0.003 

-0.068 

-0.068 

-0.003 

-0.003 


Using the steady state currents shown in Table III, one can calculate the 
instantaneous stator and rotor phase currents and torque from eqns. (29) to (33). 



ias = iqs = ^ (i qs +ke + i qs ' ke )^(k* (o*t) + (i - i^)sin(**<y e *0] 

(29) 


*ds = Z[(ids +ke + ids ~ ke )cos(k*a> e *t) - (i qs +ke - i qs - ke )sin(k*co e *t)] 

(30) 


iqr = Z[(i qr +ke + i qr ' ke )cos(k*CQ e *t) + (i^. +ke - i^r ke )sm(k*co e *t)] 

(31) 


i'dr= Z[(idr +ke + idr’ ke )cos(k* 0) e *t) - (i qr +ke - i qf T ke )sin(k*co e *t)] 

(32) 


T-M* (m/2) * (P/2) * [i qs * 4 - ids * i'j. 

(33) 

where 

P = number of poles = 2 
m = number of phases = 3, 

A/ = X m from impedance matrix 



Substitution of eqns. (29) through (32) into (33) yields a torque expression which 
contains a constant term and second harmonic components of the source constituent harmonics 
superimposed on it. The expanded torque expression appears in eqn. (34). 
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T = M (m/2) (P/2) S { ( i +ke qs *i +ke d r + i' ke qs *i' ke dr ) * ( ^ds^qr + ^ds^qr ) 
+[ ( i +ke n ^*T ke dr + i ke as*i + ^ e dr ) * ( i +ke ’ *’~ ke rke • ** +ke 


+[ ( i 


:+ke 


qs * ar ' * <75 * or' as • as * qr ■ 

y * ' ds*( ke qr +rke ds* i+ / ke qr ) 1 cos ( 2 /: ay ) 

+ i+ke ds* i ke dr ) ’ ( rke ds* i+ke dr + * ke qs* i+k£ qr ) ]sin(2/:o> e r ) (34) 


*;-ke 

qs * < 7 r 


Computer simulation results of the instantaneous phase A current and electrical torque 
are depicted in Figs. IV- 1 and IV-2. The plots in Fig. IV- 1 pertain to a comparison between the 
simulated (b) and published (a) results. The simulated plots are obtained from values in Table III 
and the expressions in eqns. (29) and (33). The plots shown in Fig. IV- 1(a) are taken directly 
from reference^]. As pointed out earlier, the calculated d-q harmonic current components match 
the published data in Table III. The d-q fundamental current components used in the simulation of 
Fig IV- 1(a) are those from reference [4]. Satisfactory correlaton is evident between the simulated 
and published plots. The torque pulsation is due to the harmonic content of the distorted current. 
The slight discrepancy in the peak values is attributed to truncation in the number of harmonic 
source voltages included in the simulation. Of particular interest is the absence of a dc offset in the 
torque plots. This is inconsistent with the torque eqn. (34). 

A repeat of the simulated results is illustrated in Fig. IV-2(b). Here, the simulated d-q 
fundamental current components, rather than the published values, are used in eqn. (29). These 
results are compared with those from reference[ 6 ], depicted in Fig. IV-2(a). The source voltage 
magnitudes from reference[ 6 ] are neither the same as those from reference [5] or those used in the 
simulation. 

It is evident that the basic form of the results from the simulation and reference[ 6 ] are 
very similar. Furthermore, the torque plots from the simulation andreference[ 6 ] clearly 
demonstrate the d.c offset embodied in the torque expression of eqn. [34]. This suggests a 
possible error in the reported d-q fundamental current components in reference[5]. 
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-IjO 

-uo 


rVvV v 


■ 0-04 SEC. 



*qs current or current (pu) 



(a) Results from Reference [4] 

Torque (pu) 


0.4 


0.2 


0 Hll III 


- 0.2 


- 0.4 


0.005 

(b) Results from Simulation 



0.015 

TIMC(S£C) 


0.025 


FIG. IV-1 - Synchronous Speed Operation Using Multiple Reference Frame 
Approach and Six-Step Inverter Induction Motor Drive 
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(a) Results from Reference [6] 


iqs current or current (pu) 



0.015 
TIME (SEC) 


2.2 


Torque (pu) 


1 .8 


1 .6 


1 .4 



0.025 0.03 


0.005 0.01 0.015 002 0.025 0.03 

TIME(SEC) 


(b) Results from Simulation 

FIG. IV-2 - Results of Six-Step Inverter Induction Motor Drive Simulation 
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IV.3. 


PWMI - Induction M achine Drive 

The simulated per-unit harmonic voltages, using the analytic relations outlined in 
Section II. 1 are summarized in Table IV for two operating frequencies, 12 and 54 Hz. The 
resulting multiple reference frame currents for the 12Hz operation are shown in Table V. 


\ Freq. W 
\(Hz) 

Harmonic \ 
Order, k \ 

12 

54 

Per-Unit 

Line-Neutral Voltages 

1 

0.221 

1.002 

5 

0.046 

0.272 

7 

0.034 

0.085 

11 

0.024 

0.091 

13 

0.022 

0.121 


TABLE IV - Computer Results of Synthesized Pulse Width Modulation Inverter 
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«(+-k) 


i ne 
1 qs 


/ne 
1 qr 

+1 

0.12860 

0.2606 

-0.0311 

-0.2676 

-5 

-0.0825 

0.0096 

0.08120 

-0.0096 

+7 

0.04380 

0.0040 

-0.0431 

-0.0040 

-11 

-0.0200 

0.0011 

0.01970 

-0.0011 

+13 

0.01570 

0.0008 

-0.0154 

-0.0008 


TABLE V - PWMI Steady-State Per-Unit Currents in Multiple Reference Frame 

(12Hz Operation) 


Substituting the steady state currents from Table V into eqns. (29) - (33) gives the simulated 
instantaneous currents, and torque in Fig. IV-4 for the 12Hz operation. The corresponding plots 
for the 54Hz operation are shown in Fig. IV-5. 

Fig. IV-4 shows the instantaneous current and torque for the 12Hz case, while fig. 
IV-5 shows the same graphs for the 54Hz case. 

Both Figs. IV-4 and IV-5 display similar harmonic distortion in the current 
waveforms, and pulsations in the torque. However, no direct comparison can be made with the 
six-step drive results, since no direct correspondence exists in the resulting voltages of the two 
drives simulated. Hence, no attempt has been made to put the plots on the same scale for 
comparison. 


current or current (pu) 



O 0.005 0.01 0.015 0.02 0.025 0.03 

TtME(SEC) 


Torque (pu) 



TIME(SEC) 


FIG. IV-4 - Steady-State Operation Using PWMI : 12Hz 
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FIG. IV-5 - Steady-State Operation Using PWMI : 54Hz 


IV-4 PDMI - Induction Ma chine Drive 

Table VI shows the results of computer simulations of the analytical voltage 
expressions derived in section II.3. The steady-state currents in the multiple reference frame are 
shown in Table VII. These results relate to a high input frequency which is approximately ten 
times the low output frequency. For this isolated PDMI operation. Table VI shows relatively high 
harmonic voltage content in the output of the inverter. This high harmonic content is reflected in 
the current and torque waveform in Fig. IV-6. This observation contradicts the normally-held view 
of relatively low harmonic content of PDM inverters. However, it must be noted that this brief 
analysis has not rigorously developed analytic relations for the PDMI in a manner similar to that of 
the PWMI. Also, the analysis has not been extended to higher output frequencies of the PDMI, 
which are more than multiples of ten of the input low frequency. 


TABLE VI - Computer Results of Pulse Density Modulated Inverter 
[Derived from waveform in reference[5]] 

[ 0) Hf - 10 ©Lf ] 


Harmonic Order 

1 

5 

7 

11 

13 


Per-Unit Voltage 

0.585 

0.153 

0.067 

0.205 

0.238 
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TABLE VII - PDMI Steady-State Per-Unit Currents in Multiple Reference Frame 

[ (Ojjf = 10 ©Lf ] 


n (+-k) 

i M ds 

ne 

1 qs 

i' ne dr 

, / he 
1 qr 

+1 

0.34060 

0.6905 

-0.0823 

-0.7091 

-5 

-0.2760 

0.0320 

0.27170 

-0.0320 

+7 

0.08670 

0.0080 

-0.0854 

-0.0080 

-11 

-0.1705 

0.0092 

0.16790 

-0.0092 

+13 

0.16690 

0.0081 

-0.1643 

-0.0081 


iqs current or current (pu) Torque (pu) 



TIME (SEC) TIME (SEC) 


FIG. IV-6 - PDM Results 
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SUMMARY 


JL- 


This report is a documentation of a computer algorithm developed to analyze 
variable-frequency A.C. drives. In particular, it contains a synthesis of a pulsewidth modulated 
inverter(PWMI), and brief descriptions of a six-step inverter and a pulse density modulated 
inverter(PDMI) drives. 

The PWMI source is generated in the algorithm. However, the six-step and PDM 
inverter sources are simulated by stating their harmonic voltage components obtained from Fourier 
analysis. A symmetrical induction machine is modeled in the arbitrary, synchronous and multiple 
reference frames. These models give the flexibility of selecting the most appropriate model for the 
intended performance prediction. 

The algorithm has been verified by a satisfactory correlation between its simulated 
results and selected data from publications. The harmonic content of the drive sources manifests 
itself in the form of distorted line currents and pulsating torque. The level of harmonic distortion is 
determined by the relative magnitudes of the source harmonic voltage magnitudes. 

In its present form, the key area of application of the algorithm is the steady-state 
performance prediction of induction machine-drives. Planned future extensions of the algorithm 
include the following: 

-synthesis and development of analytic expressions for the PDMI and analysis of 
the effects of various operating frequencies on the A.C. drive performance. 

-automatic transformation of source voltages into any of the standard reference 
frames, for subsequent hook-up to desired A.C. machine model. 

-inclusion of other performance parameters, such as operating speed profiles, and 
total harmonic distortion(THD) of machine currents. 

-inclusion of other A.C. machinery models, for possible trade-off studies towards 
space power applications. 
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